Measurement of Transverse Single-Spin Asymmetries for Di-Jet Production in 
Proton-Proton Collisions at y/s = 200 GeV 
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We report the first measurement of the opening angle distribution between pairs of jets produced 
in high-energy collisions of transversely polarized protons. The measurement probes (Sivcrs) corre- 
lations between the transverse spin orientation of a proton and the transverse momentum directions 
of its partons. With both beams polarized, the wide pseudorapidity ( — !<'?< +2) coverage for jets 
permits separation of Sivers functions for the valence and sea regions. The resulting asymmetries 
are all consistent with zero and considerably smaller than Sivers effects observed in semi-inclusive 
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deep inelastic scattering (SIDIS). We discuss theoretical attempts to reconcile the new results with 
the sizable transverse spin effects seen in SIDIS and forward hadron production in pp collisions. 



PACS numbers; 13.75.cs, 13.87.-a, 13.88. +e, 24.70.+S 

Hard scattering of light quarks has little sensitivity to 
one quark's spin orientation transverse to the scattering 
plane, due to helicity conservation (chiral symmetry) in 
the limit of zero quark mass for both quantum chromody- 
namics (QCD) and electrodynamics. Nonetheless, sizable 
sensitivity to the transverse spin of a proton has been ob- 
served at high energies in both semi-inclusive deep inelas- 
tic scattering (SIDIS) of electrons ^ and proton-proton 
collision processes with cross sections well described by 
perturbative QCD (pQCD) 0. Theoretical interpreta- 
tions of these results 3] attribute them to a combina- 
tion of soft QCD spin-dependent features of the proton 
wave function and of the final-state fragmentation of the 
struck quark into a hadron jet. Experiments that can 
unravel these contributions are essential to understand 
high-energy hadron spin dynamics. 

Of particular interest, since it arises from orbital con- 
tributions to the proton spin is the Sivers effect 
a correlation ((sp • {pp x fc^)) ^ 0) of initial-state par- 
ton transverse momentum (fcr) with the proton's spin 
(sp) and momentum (pp). This three- vector correlation 
evades time-reversal violation when orbital components 
of the proton's light-cone wave function combine with 
initial (ISI) and/or final-state interaction (FSI) contri- 
butions to the scattering process 0, 0] ■ In the spirit of 
pQCD factorization of hadron cross sections, the Sivers 
effect involves parton distribution (Sivers) functions that 
depend on both kx and longitudinal momentum fraction 
xb- In contrast to ordinary factorization, gauge invari- 
ance demands that Sivers functions incorporate pQCD- 
calculable, but process-dependent, "gauge link factors" 
describing the partonic ISI/FSI. These lead to a predicted 
sign change between SIDIS and Drell-Yan processes 

A non-zero Sivers effect revealed [s*] in SIDIS pion pro- 
duction from a transversely polarized proton target can 
be fitted with Sivers functions of opposite sign and dif- 
ferent magnitude for u vs. d quarks [3]. This account 
can be tested by treating within a common framework 
Sivers asymmetries measured for other pQCD processes, 
such as jet production in pp collisions [9|. For colliding 
proton beams moving along the ±z-axis and vertically 
(±y) polarized, the Sivers effect gives a preferential side- 
ways (±i) kinematic boost to jet momenta, causing 
a spin-dependent average deviation from 180° azimuthal 
opening angle between jets from a hard two-body parton 
scattering. We report the first measurement of this di-jet 
asymmetry, which probes gluon, as well as quark, Sivers 
functions. The data were taken in 2006 with y/s — 200 
GeV transversely polarized proton beams at the Rela- 
tivistic Heavy- Ion Collider (RHIC), providing 1.1 pb~^ 



of luminosity integrated by the STAR detector 

Continuous operation of two Siberian snakes [12| in 
each RHIC ring guaranteed that the beam polarizations 
were vertical at STAR. The spin orientation alternated 
for each successive bunch of one beam and for each pair 
of bunches of the other. Four distinct alternation pat- 
terns were used for different beam stores to minimize false 
asymmetries from accidental correlations of beam proper- 
ties with bunch number. Beam polarizations, monitored 
during each store by proton-carbon Coulomb-nuclear in- 
terference polarimeters [13], averaged 59% (57%) for the 
+z {—z) beam for this analysis, with statistical uncer- 
tainties far smaller than the ±12% relative uncertainty 
in the (online) polarimeter calibration. 

The detector subsystems critical to the present mea- 
surements are the barrel (BEMC) and endcap (EEMC) 
electromagnetic calorimeters [11], with full azimuthal 
((/)) coverage spanning pseudorapidities \r]\ < 0.98 and 
1.08 <ri< 2.0, respectively. The EMC's are subdivided 
into towers that subtend small regions in A77 and A0. 
Tower gains are calibrated, to a precision « ±5% to date, 
with minimum-ionizing particles and electrons tracked 
with star's time projection chamber (TPC). Digitized 
tower signals are summed in STAR trigger hardware over 
A77 X A(/) w 1.0 X 1.0 "jet patches". The hardware trig- 
gers used required: (1) a transverse energy sum Et > 4.0 
GeV for at least one BEMC or EEMC jet patch; (2) 
£;*?* > 14 GeV summed over the full EMC; and (3) co- 
incident signals indicating a valid collision from forward 
(3.3 < < 5.0) beam-beam counters (BBC) at each end 
of the STAR detector [13] . A software (level 2) trigger 
then passed only that subset of events with at least two 
localized (to A77 x A(j) = 0.6 x 0.6) EMC energy deposi- 
tions, with Eti(2) > 3.6 (3.3) GeV and \(f>i - (f>2\ > 60°. 

The trigger selectivity for di-jets is illustrated in Fig. 
l(a-c) by EMC information from the level 2 processor. 
The azimuthal angles 0i.2 (referred to the horizontal +x- 
axis in the STAR coordinate frame) and pseudorapidities 
111,2 (measured with respect to +z) of the two jet axes are 
obtained from i?T-weighted centroids of the EMC tower 
locations in the level 2 jet clusters. The t] values use an 
event vertex determined with coarse resolution (cr^ « 30 
cm) from the time difference between the two BBC's. 
The correlation in Fig. 1(a) is dominated by intense di- 
jet ridges centered around — 02] = 180°. 

Initial-state kx is manifested in a given event by a tilt 
of the jet axes, characterized by the deviation ]0i — 4>2\ — 
180° and the di-jet bisector angle 4>b- The Sivers analysis 
combines these features in a "signed" azimuthal opening 
angle C, chosen > 180° when cos (f>b > (implying k^ > 
0) and < 180° otherwise. STAR's left-right symmetric di- 
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FIG. 1: Distributions of events that pass the STAR di-jet 
trigger with respect to (a) di-jet azimuthal angles, (b) signed 
azimuthal opening angle (("), and (c) pseudorapidity sum, all 
using EMC energies only. The pr distribution (d) is from full 
jet reconstruction for 2% of all runs analyzed. The fit in (b) 
is from fast Monte-Carlo simulations described in the text. 



Fast MC Simulations 




FIG. 2: Fast MC simulations of di-jet Sivers asymmetries: 
(a) the C, spectra for two beam spins and (b) the resulting 
yield asymmetry vs ^ for an assumed spin-dependent centroid 
shift ±(A:r) = ±100 MeV/c. (c) The spin-dependent C cen- 
troid shift (right scale) and C, -integrated analyzing power (left 
scale) vs assumed (fcy), including (solid line) and excluding 
(dashed) cj> resolution smearing of the reconstructed jet. 



jet acceptance, reflected in the C, symmetry in Fig. 1(b), 
minimizes systematic errors in our Sivers asymmetries. 

The C distribution shape is well reproduced by "fast" 
Monte-Carlo (MC) simulations discussed below. The 
peak width (cr^ « 20°) is dominated by intrinsic kr distri- 
butions of the scattering partons, with smaller (j) resolu- 
tion contributions from the use of EMC energy alone for 
partial jet reconstruction (^jEMC-fuU _ 3.9°)^ and from 
deviations between parent parton directions and even 
fully reconstructed jet axes ((jf"ii-pa'''°" = 5.0°). These 
resolutions were determined, respectively, from the data 
themselves and from simulations utilizing the PYTHIA 
6.205 event generator [l3| and GEANT [3] modeling of 
the detector response. In the first case we compared, for 
a small sample of runs, the i^, C, and rj values (the lat- 
ter yielding g-EMC-fuU _ q qj^ determined at trigger level 
and from full jet reconstruction including offline gain cal- 
ibrations and TPC tracks. Full reconstruction, following 
the approach in but with a jet cone radius of 0.6 and 
Pt threshold of 4.0 GeV, does not greatly improve the 
net parton directional resolution ^^j^'^C-parton _ g go-j^ 

Thus, the trigger-level di-jet analysis reported here is suf- 
ficient to explore initial results and their implications for 
theoretical descriptions of the Sivers effect. 

The transverse momentum {pt) distribution from full 
jet reconstruction (Fig. 1(d)) indicates dominance of par- 
tons with XT = '^PTl\fs ~ 0.05 — 0.10. The actual xb 
range probed is broad due to the r\ coverage in Fig. 1(c). 
In a leading-order parton-parton scattering interpreta- 
tion, rii+ri2 = In (a;^^ / ) . The range 2 < (?7i -I- 7]2 ) < 3 
is then primarily sensitive to x~^^ « 0.1 — 0.4, x^^ w 
0.01 — 0.04, so that the two beams provide complemen- 



tary information on valence- and sea-dominated regions. 

Fast MC simulations in Fig. 2 illustrate some Sivers 
asymmetry measures. Two-parton scattering events were 
generated with a uniform distribution in 0i (and |02 — 
01 1 = 180°) and a pT distribution reproducing Fig. 1(d). 
Each parton was given a random initial-state kT drawn 
from a model distribution centered about zero for the y- 
component, but about ±{k'x) for the a;-component in a 
polarized proton, with the sign correlated with Sp x pp 
to simulate the Sivers effect. The sum k^^ + k^^ was 
added to the initially thrown outgoing momenta to de- 
duce boosted azimuthal angles that could then be further 
smeared with a Gaussian of g-EMC-parton _ g 30^ 

For Figs. 1(b) and 2, the model kT distribution com- 
bines a Gaussian peak with symmetric exponential tails 
enhancing larger \k^'^\, as needed to reproduce the 
roughly flat C spectrum wings. Full event reconstruction 
shows these wings to be dominated by multi-jet events, 
reflecting higher-order pQCD processes, where only the 
two jets with highest EMC energy were analyzed at level 
2. With (j> smearing included, the kT distribution fitted 
to Fig. 1(b) has an rms width ((fc^'^)2)i/2 = 1.26 GeV/c, 
consistent with the trend of earlier particle correlation 
results 17] from pp collisions. The linear relationship of 
single-spin observables to (fcf.) seen in Fig. 2(c) is rather 
insensitive to details of the kT distribution shape. 

Figures 2(a,b) show that the primary Sivers manifes- 
tation is a spin-dependent C centroid shift, leading to a 
spin up vs. down yield asymmetry antisymmetric about 
C = 180°, as predicted in [10]. We sort real data into sta- 
tistically independent ^ distributions for the four beam 
spin combinations -f-t-, etc., where the first (second) in- 
dex is the sign of the y polarization component at STAR 
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for the +z {—z) beam. To compare with predictions in- 
tegrated over the kx distributions, we extract analyzing 
powers ^^^(C > averaged over ( and 0;,, by fitting 
asymmetries measured for individual | cos(/)f,| bins: 



/P±,|cos0fc|A±^(C>7r) - [r±M-l]/[r±,icl)b)+l], (1) 

where the cross-ratios r exploit the antisymmetry in Fig. 
2(b) by treating di-jet yields Nij with spin- up and C > tt 
as equivalent to spin-down, ^ < tt. For example: 



where sums extend over — z-beam spin states j = — . 
The cross-ratio eliminates the need for independent rela- 
tive luminosities for different spin combinations and can- 
cels several potential systematic errors. {P±z) \ coS(^fc| de- 
notes beam polarization components normal to the di- 
jet bisector within each | cos 0^1 bin. The factor / — 
0.85 ± 0.07 in Eq. [T] corrects for dilution of a parton-level 
asymmetry by (f> resolution smearing (compare solid and 
dashed lines in Fig. 2(c)), with an uncertainty to allow 
for model-dependence in determining / from simulations. 
The equivalent of Eq. [2] for r^z has yields Ni^{( > tt) 
(summed over -|-z-beam spin states i) in the numerator. 
This gives An>0 when kx points preferentially leftward 
for a spin-up beam, following the Madison (18j, rather 
than the opposite Trento [l^ convention used in 0]. 

The measured asymmetries, integrated over |^ — 7r| < 
68°, are compared to calculations 20] in Fig. 3. The 
systematic error bands combine in quadrature the / un- 
certainty and the effect of multi-jet contributions to the 
C distribution wings. Limits on the latter effect are de- 
duced by looking for variations in r±z, beyond statisti- 
cal fluctuations, when we extract yields alternatively by 
changing the C integration range or subtracting a con- 
stant baseline fitted to the C, wings independently for each 
spin state. We neglect much smaller instrumental asym- 
metries from bunch-to-bunch variations in beam path or 
in azimuthally localized beam background. 

The measured asymmetries are consistent with zero, 
and remain so for higher software EMC Et thresholds. 
BBC yields analyzed with the same code reproduce the 
associated non-zero asymmetry jlJl in both magnitude 
and sign. Our results are an order of magnitude smaller 
than 7r+ SIDIS Sivers asymmetries for predominant 
di-jet sensitivity (see Fig. 3(a-b)) to both high-xs quarks 
{Ap{rii +T]2 > 2)) and low-xs gluons (A^^(r/i + m > 
2)). The ry-integrated sample (2.6 x 10^ di-jet events) has 
mean {A^^(( > tt)) = within statistical uncertainties 
~ ±0.002, probing (see Fig. 2(c)) Sivers (fcf.) preferences 
as smah as - ±3 MeV/c, or ±0.2% of ((fc^'^)^)^/^- 

Recent theory breakthroughs 2^, 2^ and our prelimi- 
nary results [2J] have stimulated rapid evolution in treat- 
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FIG. 3: Measured and calculated asymmetries vs. di-jet pseu- 
dorapidity sum for +z (left) and —z (right) beams. (a,b): 
Fraction of the calculated di-jet cross section with a quark 
(gluon) from the +z { — z) beam. (c,d): Unweighted asymme- 
tries compared with pQCD calculations [131 (histograms) for 
two models of quark Sivers functions fitted to SIDIS results 
(e,f): Asymmetries for | sin (^|- weighted yields, compared 
with calculations [20l . [2H based on twist-3 quark-gluon corre- 
lations. Vertical (horizontal) bars on the data indicate statis- 
tical uncertainties (bin widths). The systematic error bands 
exclude a ±12% beam polarization normalization uncertainty. 



ments of transverse single-spin asymmetries (SSA). Bac- 
chetta et al. [i^] deduced the gauge link structure for 
hadron or jet production in pp collisions, where both ISI 
and FSI contribute, with opposite phases. Ji et al. (23j 
demonstrated strong overlap between Sivers effects and 
twist-3 quark-gluon correlations (QGC) The pQCD 
calculations [20'| in Fig. 3 exploit these developments to 
incorporate cancellations that were absent or less severe 
in predictions made before the measurements. The 
calculations use one set of unpolarized distribution func- 
tions, yielding the parton contribution fractions in Fig. 
3(a-b), but three different models of u- and d-quark Sivers 
hmctions in Figs. 3(c-f). All assume zero gluon Sivers 
function. They are integrated over a pt range (5-10 
GeV/c) well matched to our data, and further over the 
STAR 77 acceptance [20]. We have reversed the sign of 
the calculated An to apply the Madison convention. 

The calculations in Fig. 3(c-d) use [2^ quark Sivers 
functions fitted [§] to SIDIS data Q with the c?-quark 
functional form tied either to u{xb) (VYl) or d{xB) 
(VY2) unpolarized distribution functions. For 771 ±772 ^ 2 
the A'lf^ predictions reflect the sizable HERMES asym- 
metries, diluted [iO] by partial u vs d and ISI vs FSI (the 
latter were missing in [3|) cancellations, while w 
because gluon Sivers effects are ignored. 

Figure 3(e-f) compares A^ measured and calculated 
(20! with yields in Eq. (2) weighted by JsinC] [Ilj, as 
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needed to connect to a more robustly interpretable gauge 
link structure [22] , given the apparent breakdown of fac- 
torization for back-to-back dijets [2^. The measure- 
ments, consistent with zero at all C, are hardly affected 
by the weighting, but the calculations sample a different 
Sivers function moment that can no longer be constrained 
by unweighted SIDIS asymmetries. Taking constraints 
instead from QGC fits [27|] to for inclusive forward 
hadron production in pp coUisions [1,[2^ gives di-jet An 
comparable in magnitude to our data, via more complete 
ISI vs. FSI and u vs. d cancellations |20;]. The u — d can- 
cellation can be tested in the future by filtering quark 
flavors with the leading hadron's charge sign for each jet. 

In summary, we report the first measured spin asym- 
metries for di-jet production in pp collisions. The anal- 
ysis searches for a spin-dependent sideways tilt of the 
di-jet axes sensitive to Sivers correlations between the 
proton's transverse spin and transverse momentum pref- 
erences of its partons. All measured asymmetries are con- 
sistent with zero, whether dominated by partons in the 
valence or sea regions. Perturbative QCD calculations 
can reconcile these results with sizable SSA observed for 
forward hadron production in pp and for semi-inclusive 
deep inelastic scattering via cancelling contributions from 
u and d quarks and from initial- and final-state interac- 
tions. These data constrain unified theoretical accounts 
for transverse SSA in hard pQCD processes, and their 
connection to parton orbital momentum. 
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